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Abstract 

We introduce the study of the ant colony house-hunting problem from a distributed com¬ 
puting perspective. When an ant colony’s nest becomes unsuitable due to size constraints or 
damage, the colony must relocate to a new nest. The task of identifying and evaluating the 
quality of potential new nests is distributed among all ants. The ants must additionally reach 
consensus on a final nest choice and the full colony must be transported to this single new nest. 
Our goal is to use tools and techniques from distributed computing theory in order to gain 
insight into the house-hunting process. 

We develop a formal model for the house-hunting problem inspired by the behavior of the 
Temnothorax genus of ants. We then show a fl(log n) lower bound on the time for all n ants to 
agree on one of k candidate nests. We also present two algorithms that solve the house-hunting 
problem in our model. The first algorithm solves the problem in optimal O(logn) time but 
exhibits some features not characteristic of natural ant behavior. The second algorithm runs in 
0(k log n ) time and uses an extremely simple and natural rule for each ant to decide on the new 
nest. 


1 Introduction 

Some recent work in distributed computing theory has focused on biological problems inspired by 
algorithmic tasks carried out by ant colonies; for example, [7, 9, 11, 18] study collaborative food for¬ 
aging and [3] models and proposes algorithms for task allocation within a colony. Often computer 
scientists study biologically-inspired algorithms with the aim of engineering better computing sys¬ 
tems. An alternative line of work uses tools developed to analyze and design distributed computer 
systems to better understand the behavior of biological systems [10]. In this paper, we follow the 
second approach in an attempt to better understand the house-hunting behavior of the Temnotho¬ 
rax ant. The essence of this process is a collaborative search, consensus decision, and relocation of 
the entire ant colony to a new home, with an emphasis on the decision making component. 

We first model the ant colony and its computational constraints, and formally define an algo¬ 
rithmic problem that it solves. The challenge is to develop a model that is simultaneously tractable 
to theoretical analysis and close enough to reality to give meaningful insights into ant behavior. 
For example: Is ant behavior in some sense optimal given their biological constraints? One hy¬ 
pothesis may be that behavior has been optimized through evolution. Lower bounds matching the 
performance of algorithms seen in nature can provide mathematical support for this hypothesis. 
Alternatively, results may show that behavior is far from optimal, suggesting the existence of hidden 
constraints or goals yet to be fully understood. Can we identify why certain behavioral patterns 
have developed and what environmental and biological constraints they are adaptations to? 

One motivation for studying the house-hunting process of Temnothorax ants is that it has 
received significant attention from biologists. Various aspects of the process have been studied, 
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including ants’ preferences and capabilities [15, 26], the general structure of the algorithms used by 
the ants [19, 23], variations due to colony size [13], and trade-offs between the speed and accuracy 
of the process [24], This work forms a wide basis of scientifically tested hypotheses that we can use 
for the foundation of our model and inspiration for our algorithms. 

From a distributed computing perspective, house-hunting is closely related to the fundamental 
and well-studied problem of consensus [12, 17]. This makes the problem conceptually different from 
other ant colony inspired problems studied by computer scientists. Task allocation and foraging 
are both intrinsically related to parallel optimization. The main goal is to divide work optimally 
amongst a large number of ants in a process similar to load balancing. This is commonly achieved 
using random allocation or negative feedback [1] against work that has already been completed. 
In contrast, the house-hunting problem is a decision problem in which all ants must converge to 
the same choice. Both in nature and in our proposed algorithms, this is achieved through positive 
feedback [1], by reinforcing viable nest candidates until a single choice remains. At a high level, our 
work is related to previous work on rumor spreading in biological populations [8]. 

1.1 The House-Hunting Process 

Temnothorax ants live in fragile rock crevices that are frequently destroyed. It is crucial for colony 
survival to quickly find and move to a new nest after their home is compromised. This process 
is highly distributed and involves several stages of searching for nests, assessing nest candidates, 
recruiting other ants to do the same, and finally, transporting the full colony to the new home. 

In the search phase, some ants begin searching their surroundings for possible new nests. Exper¬ 
imentally, this phase has not been studied much; it has been assumed that ants encounter candidate 
nests fairly quickly through random walking. In the assessment phase, each ant that arrives at a 
new nest evaluates it based on various criteria, e.g., whether the nest interior is dark and therefore 
likely free of holes, and whether the entrance to the nest is small enough to be easily defended. 
These criteria may have different priorities [15, 26] and, in general, it is assumed that nest assess¬ 
ments by an individual ant are not always precise or rational [25]. After some time spent assessing 
different nests, going back to the old nest and searching for new nests, an ant becomes sufficiently 
satisfied with some nest and moves on to the recruitment phase, which consists of tandem runs - 
one ant leading another ant from the old to a new nest. The recruited ant learns the candidate nest 
location and can assess the nest itself and begin performing tandem runs if the nest is acceptable. 

At this point many nest sites may have ants recruiting to them, so a decision has to be made in 
favor of one nest. The ants must solve the classic distributed computing problem of consensus. One 
strategy that ants are believed to use is a quorum threshold [22, 23] - a threshold of the number of 
ants in a candidate nest, that, when exceeded, indicates that the nest should be chosen as the new 
home. Each time an ant returns to the new nest, it evaluates (not necessarily accurately) whether 
a quorum has been reached. If so, it begins the transport phase - picking up and carrying other 
ants from the old to the new nest. These transports are generally faster than tandem runs and 
they conclude the house-hunting process by bringing the rest of the colony to the new nest. 

1.2 Main Results and Organization 

Our main results are a mathematical model of the house-hunting process, a lower bound on the 
number of rounds required by any algorithm solving the house-hunting problem in the given model, 
and two house-hunting algorithms. 

Our model (Section 2) is based on a synchronous model of execution with n probabilistic ants 
and communication limited to one ant leading another ant (tandem run or transport), chosen 


2 


randomly from the ants at the home nest, to a candidate nest. Ants can also search for new nests 
by choosing randomly among all k candidate nests. We do not model the time for an ant to find a 
nest or to lead a tandem run; each of these actions are assumed to take one round. 

Our lower bound (Section 3) states that, under this model, no algorithm can solve the house¬ 
hunting problem in time sub-logarithmic in the number of ants. The main proof idea is that, 
in any step of an algorithm’s execution, with constant probability, an ant that does not know of 
the location of the eventually-chosen nest remains uninformed. Therefore, with high probability, 
fl(logn) rounds are required to inform all n ants. This technique closely resembles lower bounds 
for rumor spreading in a complete graph, where the rumor is the location of the chosen nest [16]. 

Our first algorithm (Section 4) solves the house-hunting problem in asymptotically optimal 
time. The main idea is a typical example of positive feedback: each ant leads tandem runs to 
some suitable nest as long as the population of ants at that nest keeps increasing; once the ants 
at a candidate nest notice a decrease in the population, they give up and wait to be recruited to 
another nest. With high probability, within O(logn) rounds, this process converges to all n ants 
committing to a single winning nest. Unfortunately, this algorithm relies heavily on a synchronous 
execution and on the ability to precisely count nest populations, suggesting that the algorithm is 
susceptible to perturbations of our model and most likely does not match real ant behavior. 

The goal of our second algorithm (Section 5) is to be more natural and resilient to perturbations 
of the environmental parameters and ant capabilities. The algorithm uses a simple positive-feedback 
mechanism: in each round, an ant that has located a candidate nest recruits other ants to the nest 
with probability proportional to its current population. We show that, with high probability, this 
process converges to all n ants being committed to one of the k candidate nests within O(ATogn) 
rounds. While this algorithm is not optimal, except when k is assumed to be constant, it exhibits 
a much more natural process of converging to a single nest. In Section 6, we discuss in more detail 
possible modifications to the algorithm and various perturbations and faults to which it is resilient. 
Such robustness criteria are necessary in nature and generally desirable for distributed algorithms. 

2 Model 

Here we present a simple model of Temnothorax ants behavior that is tractable to rigorous analysis, 
yet rich enough to provide a starting point for understanding real ant behavior. 

The environment consists of a home nest, denoted no, along with k candidate new nests, identi¬ 
fied as rii for i € {1, • • • , k}. Each nest rii is assigned a quality q(i) € Q, for some set Q. Throughout 
this paper we let Q = {0,1}, with quality 0 indicating an unsuitable nest, and 1 a suitable one. 
Additionally, we assume that there is always at least one nest with q(i) = 1. 

The colony consists of n identical probabilistic finite state machines, representing the ants. We 
assume n is significantly larger than k , with k = 0{n/ log n). Additionally, ants know the value of 
n but not k, so the state machines may be parameterized by n but must be uniform for all k. This 
assumption is based on evidence that real Temnothorax ants and other species are able to estimate 
colony size and change their behavior in response [6, 4]. 

The general behavior of the state machines is unrestricted but their interactions with the envi¬ 
ronment and with other ants are limited to the high-level functions searchQ, go(), and recruitQ, 
defined below. Ants execute synchronously in numbered rounds starting with round 1. In each 
round, each ant may perform unlimited local computation (transition through an arbitrary sequence 
of states), along with exactly one call to one of the functions: searchQ, goQ, or recruitQ. 

At the end of each round r, each ant a is located at a nest, denoted by £(a, r) € {0,1, • • ■ , fc}; 
initially, before round 1, all ants are located at the home nest. The value of £(a,r) is set by the 
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calls to searchQ, go(), or recruit() made by the ant in round r. Let c(i,r) = \{a\£(a,r) = i}| 
denote the number of ants located in nest n* at the end of round r. 

In each round, each ant a performs exactly one call to the following functions: 

• searchQ: Returns a triple < i,q(i),c(i,r ) > where i is chosen uniformly at random from 
{1, - - - ,k}. Sets £(a,r) := i. This function represents ant a searching randomly for a nest; 
the return value of the function includes the nest index, the nest’s quality and the number of 
ants at the nest. 

• go(i): Takes input i E {1, • • • , k} such that there exists a round r' < r in which £(a, r') = i. 
Returns c(i,r). Sets £(a,r) := i. The function represents ant a revisiting a candidate nest np, 
the function returns the number of ants at nest rii at the end of round r. 

• recruit(b, i): Takes input 6 E {0,1} and i E {1, • • • ,k} such that there exists a round r' < r 
in which £(a,r') = i. Returns a pair < j,c(0,r) > where j E {1, ••• ,k}. Sets £(a,r) := 0. 
The return value j is determined as follows. Let R be the set of all ants that call recruit( •, •), 
and let P be a uniform random permutation of all ants in R. Let S C R be the set of ants 
that call recruit{ 1 , •). 


Algorithm 1: Return value j of recruit (•, •) for all ants a E I! 
l M: a set of pairs (a, aQ of ants, initially 0 7 for i = 1 to \P\ do 


2 for i 

= 1 to P do 8 

if ( a*,apu \) E M then 

3 

if ap(j) E S and (• ,a P 0 M then 9 


return j to ant a P ^ where j is 

4 


a' := uniform random ant from R io 


input to recruitil, j) called by a* 

5 


if (a', •) fL M and (-, a 1 ) 0 M then X1 

else return j to ant ap^\ where j is 

6 


| M := M U (ap(j),aQ 12 

input to recruit(-,j) called by a P ^\ 


In short, all actively recruiting ants in S randomly choose an ant to recruit. P simply serves as 
tie breaker to avoid conflicts between recruitments. It is important to note that this process is not 
a distributed algorithm executed by the ants, but just a modeling tool to formalize the idea of ants 
recruiting other ants randomly without introducing dependency chains between the ordered pairs 
of recruiting and recruited ants. Algorithm 1 can be thought of as a centralized process run by the 
environment of the home nest in order to pair ants appropriately. We believe our results also hold 
under other natural models for randomly pairing ants. 

An ant recruits successfully if it is the recruiting ant (first element) in one of the pairs in M . If 
an ant recruits successfully or is not recruited, recruit () just returns the input nest id i. Otherwise, 
it returns the id of the nest that the ant is recruited to. 

We intentionally define recruitQ to locate ants back to the home nest, where recruitment 
happens. Additionally, go(i) is only applicable to a candidate nest rti for i 7 ^ 0 , so calling recruitQ 
is the only way to return to the home nest. Since ants are required to call one of the three functions 
in each round, if an ant does not search (call searchQ) or go to a nest (call goQ), it is required to 
stay at the home nest. In this way, all ants located in the home nest participate in recruitment in 
each round, either actively recruiting (6 = 1 ) or waiting to be recruited (6 = 0 ). 

Recruitment encompasses both the tandem runs and direct transport behavior observed in 
Temnothorax ants. Since direct transport is only about three times faster than tandem walking 
[ 21 ], and since we focus on asymptotic behavior, we do not model this action separately. 

Next, we prove a general statement about the recruitment process that will be used in the proofs 
of our algorithms. 
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Lemma 2.1. Let a be an arbitrary ant that executes recruit( 1,-) in some round r, and assume 
c(0,r) > 2. Then, with probability at least 1/16, (a,-) € M; that is, ant a succeeds in recruiting 
another ant. 

Proof. One requirement for ant a to succeed recruiting is that it does not get recruited by another 
ant positioned earlier in the random permutation P than a. With probability 1/2 ant a is located 
in the first half of the random permutation, so conditioning on that, the probability that a is not 
recruited by some ant with higher precedence in some round r is at least: 

/ ^ \ c (0,O/2 ^ 

V 1 “ c(0, r )) - 4(V 2 ) - 2’ 

where c( 0 , r) > 2 is the number of ants at the home nest in round r. 

Another requirement for ant a to succeed recruiting is to choose to recruit an ant that has 
not itself successfully recruited and has not been recruited by the preceding ants. Recall that we 
conditioned on ant a being located in the first half of the permutation P. With probability 1/2, ant 
a chooses to recruit an ant a' located in the second half of the permutation; therefore, a' appears 
after a in the permutation and so a' has not recruited yet when a recruits. With probability at 

/ \c(0,r)/2 

least (1 — > 1 / 2 , ant a' is not recruited by any ant other than a in round r. 

Therefore, in total, the probability that ant a recruits successfully is at least: 1/2 to be located 
in the first half of the permutation, 1/2 to not be recruited by another ant, 1/2 to recruit an ant 
in the second half of the permutation, and 1 /2 for this chosen ant to not be recruited by any other 
ant. This gives a total probability of at least 1/16. □ 

Problem Statement: An algorithm A solves the Househunting problem with k nests in T £ N 
rounds with probability 1 — 6, for 0 < 6 < 1 , if with probability 1 — 5, taken over all executions 
of A, there exists a nest i € {1, • • • , k} such that q(i) = 1 and l(a, r) = i for all ants a and for all 
rounds r > T. 

3 Lower Bound 

In this section, we present a lower bound on the number of rounds required for an algorithm to solve 
the house-hunting problem. The key idea of the proof is similar to the lower bounds on spreading a 
rumor in a complete graph [16] where neighbors contact each other randomly. Assuming a house¬ 
hunting process with a single good nest, its location represents the rumor to be spread among all 
ants and communication between random neighbors is analogous to the recruiting process. 

Assume only a single nest, n w , has quality 1, and so it is the only option for the ants to relocate 
to. Additionally, assume that each ant is able to recognize nest n w once it knows its id. These 
assumptions only restrict the environment and add to the ants’ power, so a lower bound under 
these assumptions is sufficient. Note that each ant has only two ways of reaching a nest: (1) by 
searching, or (2) by getting recruited to it. Throughout this section, an ant is considered to be 
informed if it know the id of the winning nest n w ; otherwise, an ant is considered to be ignorant. 

Lemma 3.1. With probability at least 1/4 an ant that is ignorant at the beginning of round r 
remains ignorant at the end of round r. 

Proof. In round r, ant a may either (1) not participate in recruitment or search (2) participate in 
recruitment at the home nest or (3) search. In the first case, ant a does not call either recruit() 
or search() and is guaranteed to remain ignorant. 
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In the second case, ant a calls recruit(0, •). Let X r be the number of informed ants at the 
beginning of round r that are recruiting at the home nest to the winning nest; these ants are calling 
recruit(1, w). Then, the probability that ant a remains ignorant after round r is at least: 

1 \ x r ( 1 \ \ 

c(0, r) / — \ c(0, r) / ~ 4 

where c(0, r) > 2 is the number of ants at the home nest at the beginning of round r. Note that if 
c(0, r) < 2, ant a is forced to recruit itself, so it remains ignorant. 

In the third case, for k > 2, the probability that ant a remains ignorant after one round of 
searching is (k — l)/k > 1/2. 

Therefore, the overall probability that an ignorant ant remains ignorant in round r, after either 
searching or participating in recruitment, is at least 1/4 (since 1/2 > 1/4). □ 

Let random variable X r denote the number of ignorant ants at the beginning of round r. In 
order for the algorithm to solve the house hunting problem, it is necessary that X r = 0 for some 
round r. Otherwise, there is at least one ant that is not committed to the winning nest n w . 

Theorem 3.2. For any constant c > 0, let A be an algorithm that solves the HOUSEHUNTING 
problem with k > 2 nests in T rounds with probability at least l/n c . Then, T = fl(logn). 

Proof. By Lemma 3.1, with probability at least 1/4 an ant that is ignorant at the beginning of 
round r remains ignorant at the end of round r. The probability that an ignorant ant remains 
ignorant for r = (log 4 n)/2 — log 4 (12c) consecutive rounds is at least (l/4) r > 12 c/y/n, where c > 1 
is an arbitrary constant 1 . So the expected number of ignorant ants at the beginning of round 
r = (log 4 n)/2 — log 4 (12c) is E[X r ] > (12 c/y/n)(n) = Ylcyfn. 

Note that the random variable X r can be expressed as a sum of indicator random variables 
Xf. where for each ant a, X°: is 1 if ant a is ignorant, and 0 if it is informed, at the beginning 
of round r. Unfortunately, these indicator random variables are not independent, so we cannot 
directly apply a Chernoff bound. However, we can still show that their sum is bounded by defining 
random variables that dominate the Xf variables. 

Let random variable Yf be a random variable such that P[Yf = 1] = (l/4) r for each a and 
each r. Let Y r be the sum of these independent random variables for all ants a. Note that 
E[Y r ] = n(l/4) r , so for r = (log 4 n)/2 —log 4 (12c), E[y r ] = Ylcyfn. Therefore, by a Chernoff bound: 

P[Y r < 6cy/n\ < e -( 12c ^)/ 12 < e~ clnn < —. 

n c 

By Lemma 3.1, the probability that an uninformed ant remains uninformed is at least 1/4, 
regardless of the other ants’ actions in round r. Therefore, P[Xf = 11 all Xf for a' before a in P] > 
(l/4) r , where P is the random permutation used to model the recruitment process. So, by the 
definition of Yf : 

P[Xf = 11all Xf for a' before a in P] > Q = P[Y r a = 1] 

By Lemma 1.18 in [5][Chapter 1] (also stated as Lemma A.l in the Appendix), P[X r < x] < 
P[Y r < x\ for any x < n. In particular, for r = (log 4 n)/2 — log 4 (12c), it follows that : 

P\X r < 6cy/n\ < P[Y r < Gcy/n\ < — 

n c 

1 We assume a constant is asymptotically smaller than y/n. 
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Therefore, with probability at least 1 — l/n c , at least 6 Cy/n ants are ignorant at the end of 
round r = (log 4 n)/2 — log 4 (12c). These ants are not informed of the id of the winning nest, and 
so cannot be located at this nest. 

Since algorithm A solves the HOUSEHUNTING problem with probability at least 1 /n c in T 
rounds, then, with probability at least l/n c after T rounds, all ants are informed of the winning 
nest. We showed that with probability at least 1 — l/n c , at least 6 Cyfn ants are ignorant at the end 
of round r = (log 4 n)/2 — log 4 (12c). Therefore, T > r = 12 (log n). □ 

4 Optimal Algorithm 

We present an algorithm that solves the Househunting problem and is asymptotically optimal. 
In the key step of the algorithm, each ant tries to recruit other ants to the nest it found after 
searching; after each round of recruiting, each ant checks if the number of ants at its nest has 
increased or decreased. Nests with a non-decreasing population continue competing while nests 
with a decreasing population drop out. In each round, the population of at least one nest is non¬ 
decreasing, so at least one nest will remain in the competition. Additionally, other nests drop out 
at a constant rate, meaning a single winning nest will be identified quickly. 

This algorithm relies heavily on the synchrony in the execution and the precise counting of the 
number of ants at a given nest, which makes it sensitive to perturbations of these values, and there¬ 
fore, not a natural algorithm that resembles ant behavior. However, the algorithm demonstrates 
that the HouseHunting problem is solvable in optimal time in the model of Section 2. 

4.1 Algorithm Pseudocode and Description 

The pseudocode of the algorithm is presented in Algorithm 2. Each call to the functions from 
Section 2 (in bold) takes exactly one round. The remaining lines of the algorithm are considered 
to be local computation and are executed in the same round as the preceding function call. 

Throughout the algorithm’s execution, each ant is in one of four states: search, active , passive, 
or final, initially search. Based on the state of the ant, it executes the corresponding case block 
from the pseudocode. An ant is said to be committed to a nest n* if nest = i. 

Note that the search subroutine is executed only once, during the first round, and the final 
subroutine is what active ants do at the end of the execution to recruit all ants to the winning nest. 
The other two subroutines, active and passive, represent the actions of active (recruiting) ants, 
and ants from bad/dropped-out nests, respectively. Each one of these subroutines takes exactly 
four rounds; each round is labeled as either Rl, R2, R3, or R4 in the pseudocode. The subroutines 
are carefully scheduled in such a way that these two types of ants do not meet until the end of 
the competition process when a single winning nest remains; otherwise, the competition between 
ants from competing nests would be slowed down by ants from dropped-out nests. Therefore, the 
active and passive subroutines are padded with recruit(0 , •) and go{nest ) calls to achieve such 
interleaving (lines 13, 18-19, 35-36, 42 are such padding rounds). 

• Search (lines 6—11): In the first round, each ant searches for a nest and commits to it. 
If the nest has quality of 0, the ant moves to the passive state; otherwise, it moves to the 
active state. 

• Passive (lines 12—19): A passive ant is an ant committed to either a bad or a dropped-out 
nest. The ant spends a round at its (non-competing) nest, then it goes home to be recruited. 
This call to recruit(0,nest) never coincides with a recruit[\,nest) of an active ant, so a 
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Algorithm 2: Optimal HouseHunting Algorithm 


1 state: {search, active, passive, final}, 

2 initially search 

3 nest : nest index i E {0, • • • , k}, initially 0 

4 count : an integer in {0, • • • , n}, initially 0 

5 quality : a boolean in {0,1}, initially 0 


6 

7 

8 
9 

10 

11 


case state = search 

< nest, count, quality >:= searchQ (Rl) 
if quality = 0 then 
| state := passive 

else 

I state := active 


12 

13 

14 

15 

16 

17 

18 
19 


case state = passive 
go(nest) (Rl) 

< nest t ,- >:= recruit(0, nest) (R2) 
if nestt ~f~ nest then 
nest := nestt 
state := final 
go(nest) (R3) 
go(nest) (R4) 


20 case state = final 

21 I < nest, ■ >:= recruit(1, nest) (Rl) 


22 case state = active 


23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 


< nestt, - >■= recruit(1, nest) (Rl) 
countt := go(nest t ) (R2) 
if ( nestt = nest) and 
(countt > count ) then 
count := countt 
go(nest) (R3) 

< -,counth >:= recruit(0, nest) (R4) 
if counth = count then 
| state := final 
else if ( nestt = nest) and 

(i countt < count ) then 
state := passive 
recruit(0, nest) (R3) 
go(nest) (R4) 
else 

nest := nestt 
count n := go(nest) (R3) 
if count n < countt then 
| state := passive 
go(nest) (R4) 


passive ant can only get recruited by an ant in the final state calling recruit(\,nest). Once 
successfully recruited, the passive ant moves to the final state and commits to the new nest. 

• Final (lines 20—21): An ant in the final state is aware that a single winning nest remains, 
so it recruits to it in each round. This call to recruit(l,nest ) coincides with the call to 
recruit}0 , nest ) of passive ants, so once a single nest remains, passive ants are recruited to it 
in every fourth round. 

• Active (lines 22—42): An active ant tries to recruit other ants to its competing nest by 
executing recruit(l,nest). The ant then goes to the resulting nest to count the number of 
ants there. Based on the values of the resulting nest (nestt) and count (countt), we consider 
three cases: 

— Case 1 (lines 25-31): nestt is the same as the nest that the ant is committed to and the 
number of ants in that nest has not decreased; therefore, the nest remains competing. 
As a result, the ant updates the new count and spends an extra round at the nest that 
has a special purpose with respect to Cases 2 and 3 below. Finally, the ant checks if the 
number of ants at the home nest is the same as the number of ants at the committed 
nest; if this is the case, it means that all ants have been recruited to a single winning 
nest and the ant switches to the final state. 

— Case 2 (lines 32-36): nestt is the same as the nest that the ant is committed to but the 
number of ants has decreased; therefore the nest is about to drop out. As a result, the 












ant sets its state to passive and spends a round at the home nest, which coincides with 
the round an active ant spends at the committed nest in Case 1. 

— Case 3 (lines 37-f2): nestt is different from the nest that the ant is committed, which 
indicates the ant got recruited to another nest. Although it already knows the number of 
ants ( countt ) at the new nest, the ant updates that count ( count n ). The reason for this is 
to determine whether this new nest is about to compete or drop out. If countt = count n , 
the nest is competing because the active ants in Case 1 are spending the same round 
at the committed nest; if countt. > count n , the nest is dropping out because the ants in 
Case 2 already determined a decrease in the number of ants and are spending this round 
at the home nest. 

4.2 Correctness Proof and Time Bound 

As written, Algorithm 2 never terminates; after all ants are committed to the same nest and in the 
final state, they continue to recruit (each other) in every round. This issue can easily be handled 
if ants check whether the number of ants at the home nest is the same as the number of ants at the 
candidate nest. However, for simplicity, we choose not to complicate the pseudocode and consider 
the algorithm to terminate once all ants have reached the final state and, thus, committed to the 
same unique nest. 

Proof Overview: The correctness proof and time bound of Algorithm 2 are structured as 

follows. By defining a slightly different but equivalent recruitment process, we show, in Lemma 
4.1, that a competing nest is equally likely to continue competing and to drop out. Consequently, 
as we show in Lemma 4.2, each competing nest has a constant probability of dropping out. We 
put these lemmas together in Theorem 4.3 to show that, with high probability, Algorithm 2 solves 
the HOUSEHUNTING problem in O(logn) rounds: O(logfc) rounds to converge to a single nest and 
O(logn) rounds until all passive ants are recruited to it. 

Let be the set of round numbers r such that r e R% iff r mod 4 = 3; these rounds are 
exactly the rounds in which only ants committed to active nests are located at the home nest and 
try to recruit each other. Also, let C{i,r) for each r € R 3 and each nest n % denote the set of ants 
committed to nest rq and competing at the home nest. This implies that | C(i,r)\ = c(0,r) 

for r £ R 3 . 

Let random variable Xf, for each ant a and each round r £ R 3 , take on values —1, 0 or 1 as 
follows. If ant a gets recruited by another ant in round r, then Xf = —1; if ant a successfully 
recruits another ant, then Xf = 1; otherwise, Xf = 0. 

Let random variable Yf denote the change in the number of ants at nest n* after one round of 
recruiting, where r £ R 3 : 

\ri \ ' va 

1 r / / ■ 

a£C(i,r) 

Informally speaking, the change in population of nest nt is simply the sum of identically dis¬ 
tributed {—1,0,1} random variables with mean 0 that take on non-zero values with constant proba¬ 
bility. Therefore, the sum of these variables is negative with constant probability. However, proving 
this fact requires a more rigorous argument because the Xf variables are not independent. 

As specified by the recruitment model, Yf is the result of a randomized recruitment process 
in which all ants are ordered by a random permutation. Then, in order of the permutation, ants 
choose uniformly at random other ants at the home nest to recruit. The random variables involved 
in this process are the random permutation P as well as the set of random choices of ants. Consider 
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the vector o r , of length c( 0 ,r), where o r (a) = — 1 if ant a is chosen by another ant, o r (a) = 1 if 
ant a chooses an ant a' such that a' had not already chosen an ant and a' was not already chosen 
by another ant, and o r (a) = 0 otherwise. 

Lemma 4.1. Let nest n* be a competing nest in round r E P 3 . Then, P[Y, / < 0] = P[Yf > 0]. 

Proof. Random variable Yf is simply the sum of values of o r (a) corresponding to ants a € C(i,r). 
Since o r has exactly the same number of —Is and Is, it is possible to choose some (non-random) 
permutation P 1 that swaps the locations of the — Is and Is in o r . Choosing a random permutation P 
in the recruitment process and then applying the swapping permutation P 1 negates the value of Yf. 
Moreover, choosing a uniform random permutation P and then applying this swapping permutation 
P' still yields a uniform random permutation. So, —Yf is distributed according to the exact same 
distribution as Yf. Therefore, Yf is symmetric around 0 and P[Yf < 0] = P[Yf >0]. □ 

Lemma 4.2. Let nest Hi be a competing nest in round r E R 3 . If \C(i,r)\ < c(0,r), then P[Yf < 

0 ] > 1 / 66 . 

Proof. Let Pi and P 2 be uniform random permutations. Consider using Pi as the uniform permu¬ 
tation in in the recruitment process that determines Yf. Permutation P 2 swaps the position of a 
fixed ant a* not committed to nest Hi with the position of an ant a chosen uniformly at random 
among the ants C(i,r) (the ants committed to nf). The existence of such an ant is guaranteed by 
the assumption that \C(i,r)\ < c(0, r). 

By Lemma 2.1, P[Xf = 1] > 1/16. Conditioning on Yf = 0, at most 1/2 the ants a E C{i,r) 
have Xf = 1. Therefore, with probability at least 1/2, the ant a chosen uniformly at random by 
permutation P 2 has Xf < 1. Conditioning on Xf* = 1 can only increase this probability. So with 
probability at least (1/2) (1/16), applying the composition of Pi and P 2 to compute Yf increases 
its value by at least 1 . Since the composition of Pi and P 2 is also a uniform random permutation, 
the distribution of Yf remains exactly the same as the case when only Pi is applied. Therefore, 

p[y; = 0] < 1 - y 2 P[Yf = 0] 

99 

PK' = 0] < 33. 

By Lemma 4.1, P[Yf < 0] = P[Yf > 0]. Therefore, P[Yf < 0] > 1/66. □ 

Theorem 4.3. For any constant c > 0, with probability at least 1 — 1 /n c , Algorithm 2 solves the 
Househunting problem in O(logn) rounds. 

Proof. In the first round, all ants search for nests, so the expected number of ants located each 
good nest is n/k. Assuming k < n/ (12(c + 1) log n) = 0{n/ log n), by a Chernoff bound, it follows 
that, with probability at least 1 — l/n c+1 , at least n/( 2 k) ants visits each good nest. 

Let k r be a random variable denoting the number of competing nests in round r E P 3 . Suppose 
k r > 1, and so \C(i,r)\ < c(0,r). By Lemma 4.2, P[Yf < 0] > 1/66 for each nest n* among the 
k r competing nests. Therefore, E[fc r + 4 ] < (65/66)fc r . Also, note that for any round r, k r < k. 
For r = log 66 / 65 k + (c + 1) log 66/ / 65 n = O(logn), where c > 0 is an arbitrary constant, it follows 
that E[Ay] = l/ra c+1 . By a Markov bound, P[k r > 1] < l/n c+1 , so with probability at least 
1 — l/n c+1 , k r < 1. Note that in any given round, it is not possible for all nests to experience 
decrease in population, so for each round r, we have k r > 1. So, overall, after O(logn) rounds, 
with probability 1 — l/n c+1 , there is exactly one competing nest. 
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After there is exactly one competing nest in some round r £ i? 3 , all ants committed to it 
switch to the final state and start recruiting the passive ants during every fourth round. Each 
recruited ant also transitions to the final state, resulting in at most O(logn) rounds until all ants 
are committed to the winning nest. 

Therefore, in total, with probability at least 1 — l/n c each nest is discovered by at least one ant 
and all ants are committed to a single good nest in O(logn) rounds. □ 

5 Simple Algorithm 

We now give a very simple algorithm (Algorithm 3) that solves the Househunting problem in 
0(k\ogn) rounds with high probability. The main idea of the algorithm is that all ants initially 
search for nests and those that find good nests simply continuously recruit to their nests with 
probability proportional to nest population in each round. Ants in larger nests recruit at higher 
rates, and eventually their populations swamp the populations of smaller nests. This process is 
similar to the well-known Polya's urn model [2]. 

5.1 Algorithm Description and Pseudocode 

In each round of Algorithm 3, each ant can be in one of two states: active or passve, initially 
starting in the active state. In the first round of the algorithm, all ants search for nests; the ants 
that find good nests remain in the active state, and the ants that find bad nests switch to the 
passive state. Then, the algorithm proceeds in alternating rounds of recruitment by all ants at 
the home nest (either active recruit(1, •) or passive recruit(0, •)), and population assessment at 
candidate nests. In each round of population assessment, each ant chooses to recruit actively in 
the next round with probability count/n, where count is the assessed population at the candidate 
nest, and n is the total number of ants. When a passive ant gets recruited to a nest, it commits 
to the new nest and becomes active again. When an active ant gets recruited to a different nest, it 
commits to the new nest and starts recruiting for that new nest. 


Algorithm 3: Simple House-Hunting 

1 state : {active,passive}, initially active 

9 

case state = passive 

2 < nest, count, quality >:= searchQ 

10 


nestt : = recruit(0, nest) 

3 if quality = 0 then 

11 


if nestt 7^ nest then 

4 

state := passive 

12 



state := active 

5 case state = active 

13 



nest := nestt 

6 

b := 1 with probability 0 otherwise 

14 



count := go(nest) 

7 

nest := recruit(b, nest) 
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count := go(nest) 






5.2 Correctness Proof and Time Bound 

For each nest ni and each round r let random variable p(i,r) = c(i,r)/n denote the proportion 
of ants at nest n* in round r. By the pseudocode, ants located at nest n t in round r will recruit 
with probability p(i,r) in round r + 1. Define S 2 (r) = Yli=iP(h r ) 2 i the expected proportion of 
ants that will recruit in total, in round r. Since Yli=iP(^ r ) = C by the l\ versus I 2 norm bound 
S 2 (r) > l/k. 
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Proof Overview: The correctness proof and time bound of the algorithm are structured as 

follows. First, in Lemmas 5.1, 5.2, and 5.3, we show some basic bounds on the expected number of 
ants recruited in each round and the change of the value of p(i, r) for a single nest in each round of 
recruiting. Then, in Lemmas 5.5, 5.6, and 5.7, we use a Taylor series expansion of the ratio between 
the populations of two nests to show that in expectation this ratio increases multiplicatively by 
(l+n(l/&;)) in each round of recruiting, provided that both nests have a £1(1/k) fraction of the total 
population. On the other hand, if a nest has less than a £1(1/k) fraction of the total population, in 
Lemmas 5.8 and 5.9, we show that the ants in such a nest recruit so slowly that the nest completely 
drops out within O(klogn) rounds with high probability. Finally, in Theorem 5.11, we consider all 
( 2 ) pairs of nests to show that only a single nest will remain in each pair within O(klogn) rounds, 
ensuring, by a union bound, that a single nest remains overall. 

Throughout this section, let c > 0 be an arbitrary constant and let d be an arbitrary constant 
such that d > 64. For the sake of analysis, assume that k < yjn/ (8 d 2 (c + 6) log n) = (D(^Jn/ logn). 
While we feel that this assumption is reasonable, we are also hopeful that it could be removed. Let 
R\ be the set of all odd rounds, excluding round 1; by the pseudocode, in these rounds, ants are 
located at the candidate nests. Therefore, for each r € R\ and each ant a, £(a,r) = i ^ 0; for each 
round r' £L R\ and for each ant a, £(a, r') = 0. 

5.2.1 Change in Population of a Single Nest in One Round 

We first study how the population of a single nest changes in a single round. Intuitively, we expect a 
p(i , r) proportion of the ants in m to recruit, and a X 2 (r) proportion to be recruited. Therefore, we 
expect the population to change by p(i,r ) [p(i,r) — £ 2 (r)] • Qualitatively, we show this in Lemma 
5.3, modulo constant factors; the main technical difficulty is handling dependencies in the random 
recruiting process. Before we are ready to prove Lemma 5.3, we show two lemmas that state the 
expected outcome of a single ant recruiting. 

As before, define random variable Xf to take on value —1 if ant a is recruited away from its 
current nest in round r, 1 if it successfully recruits another ant, and 0 otherwise. 

Note that the following lemma is stated in a slightly generalized way: it applies to an ant a 
recruiting with a fixed probability p in a given round r. By the pseudocode of the algorithm it is 
clear that whenever an ant recruits, it is always with probability p(i, r), so the majority of the time 
we will apply the lemma for p = p(i,r ); however, the general statement of the lemma helps reason 
about the expected value of X/ with respect to two ants from two different nests in the proof of 
Lemma 5.2. 

Lemma 5.1. Let rii be any nest and let ant a be located in nest m in some round r € R\. Also, 
suppose ant a recruits with probability p in round r. Then, there exist functions £1 and £2 such that 
E[X“] = p£i (i, r) — £2 (i,r) and £1 (i,r) > £ for a fixed constant £ > 0. 

Proof. Let be an indicator variable indicating whether some ant a chooses to recruit (executes 
recruit(1, •)) in round r + 1. By assumption, Pr \A a r = 1] = p. Whether ant a is actually successful 
in recruiting another ant depends on: (1) the order of the random recruiting permutation P, (2) 
the choices of other ants to recruit or not, and (3) the choices by recruiting ants of whom to recruit. 
Let B be a random variable encompassing all these random variables that affect Xf, excluding 
Therefore, B is a triple ( P , {A£ |o' / a}, {1, ■ ■ ■ , n} n ) and it takes on values from the set B of all 
such triples. The expected value of Xf is: 

E[X“] = p r [B = b\ ■ (pE [X?\B = b, A a r = 1] + (1 - p)E [X/\B = b,A a r = 0]). 
fees 
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Fix some value B = b. We consider several cases based on the nine possible pairs of values of A “ 
and (A“ has two possible values and X" has three possible values). However, note that some 
permutations of these values are not allowed; that is, if = 0 (the ant chooses not to recruit), 
it is not possible that A“ = 1 (the ant succeeds in recruiting). Also, since B = b is already fixed, 
ant a is either chosen by another ant or not, regardless of the value of A“. This fact rules out two 
more cases: (1) the case where = 0 if A“ = 0, and X“ = — 1 if A“ = 1, and (2) the case where 
= —1 if = 0, and X“ = 0 if A“ = 1. In (1), since the choices of the other ants are already 
fixed and included in B = 6, it is not possible that ant a gets recruited if it chooses to recruit but it 
does not get recruited otherwise. Similarly, in (2), it is not possible that ant a gets recruited when 
it chooses not to recruit but does not get recruited when it fails to recruit. 

All remaining cases are listed below: 


Case 1: = —1 for both A“ = 0 and A“ = 1. That is, the ant is recruited by another ant no 

matter its decision. 

Case 2: = — 1 if A° r = 0, and A“ = 1 if A“ = 1. That is, if the ant chooses to recruit it succeeds 

and if not, it is recruited by another ant. 

Case 3: = 0 for both = 0 and A“ = 1. That is, whether the ant chooses to recruit or not, 

it will not be part of a successful recruitment. 

Case 4: = 0 if A“ = 0, and X.“ = 1 if = 1. That is, if the ant chooses to recruit it succeeds, 

but if not, it is not recruited by another ant. 


Let p° c be the probability of case c occurring given B = b. The expected value of X“ is: 

Pr = ' (p ( p 2 + p\ ~ Pi) + (1 - P) (-Pi - P 2 ) ) 


b&B 

= [B = b] (p (2 P 2 +P 4 ) - (p\ +p b i )) 

b&B 

= pE B \^p\ +P 4 I - 


p\+p\ 


Now, we need to show that E# \2p\ + p\] > ^ for some fixed constant £ > 0. By Lemma 2.1, 
each recruiting ant succeeds with some probability lower bounded by a constant. A recruiting ant 
only succeeds in cases 2 and 4 so we have E b\p 2 + p\] > £ for some constant £ > 0. Therefore, 
Eb[2?4 +P 4 ] > E b\P 2 +P 4 ] > 6 

Overall, E[X“] = p£i(i,r) — £ 2 (*, 7 ™), where £i(i,r) = E^[2p2 +P 4 ] is any constant greater than 
£ > 0 and £ 2 ( 2 , r) = E^[p5 +P 2 ] i s some unspecified value. □ 

Note that the values of £1 (i,r) and £ 2 ( 1 , 7 -) are slightly different for different nests because the 
behavior of an ant depends on the nest it is committed to. For example, the different probabilities 
corresponding to the cases in Lemma 5.1 vary based on the population of a given nest in a given 
round and the populations of the remaining nests in that round. Next, we prove a bound on the 
£i(i,r) and £2 (i,r) values of two different nests. 

Fix a constant £ and functions £1 and £2 such that £1 (i,r) > £ for all rounds r and all nests rij. 

Lemma 5.2. Let rii and nj be two nests with p(i,r ) < p(j,r ) in some round r € R\. Then, 
£i(z,r) -p{i,r) - £ 2 (i,r) < £1 (j,r) -p(i,r) - £ 2 (j,r). 
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Proof. Let a* and a,j be two ants located in nests n* and rij , respectively, in round r E R\. Ant a* 
recruits with probability p(i,r ) and ant aj recruits with probability p(j,r ) where p(i,r) < p(j,r). 
By Lemma 5.1, applied to ant a*, recruiting with probability p(i,r), and nest n^, it follows that 

Suppose we swap the locations of ants a* and a,j and let both ants recruit with probability 
p(i,r)', since the ants just switch locations, the populations of the two nests remain the same, so 
the rest of the ants still recruit with probability p(i , r) in nest n* and p(j , r ) in nest iij 2 . By Lemma 
5.1, applied to ant a*, recruiting with probability p(i,r), and nest rij , it follows that the expected 
value of X“ ,: in this hypothetical situation is E'fA^b] = p(i,r)£i(j,r) — £2 (j, r )- 

As a result of switching the positions of ants a* and a,j, we are essentially lowering a^'s probability 
to recruit. Therefore, the probability that ant at succeeds recruiting, if it chooses to recruit, 
increases, and the probability that ant a, is recruited, if it chooses not to recruit, decreases. 

Therefore, E[X^] < E'[X“*], and so p(i,r)£i(i,r) — £2 {i,r) < p(i,r)£i(j,r) — □ 

Next, we use Lemma 5.1 to calculate the expected change in the value of p(i, r ) after one round 
of recruiting. 

Lemma 5.3. For each nest Hi and each round r € R\: 

E \p(i, r + 2)] = p(i, r) • [1 + £1 (i, r) ■ p(i, r) - £ 2 (l r)]. 

Proof. By linearity of expectation, 

E\p(i,r + 2)]=p(i,r) + ^ £ E[A“] = p(i, r) ■ (1 + E[A“]), 

{a\e(a,r)=i} 

where in the last equation a is a fixed ant in nest n t and the equality follows from the fact that 
Xf is identically distributed for each ant in the same nest. 

By Lemma 5.1, E[X“] = p(i, r) — £ 2 (l r), so: 

E\p(i,r + 2)\ =p(i,r) ■ (1 +p(i, r)£i(i, r) -£ 2 (br)). □ 

Lemma 5.3 shows that the population change in a single nest depends quadratically on p(i,r). 
With this fact, we can show that larger nests will tend to ‘swamp’ smaller nests, causing their 
populations to drop to 0. 


5.2.2 Relative Changes in the Populations of Two Nests in One Round 

We first define a measurement of the population gap between two nests. 


Definition 1. For any two nests rii and nj, let nn(i,j,r ) E {*, j} be the id of the nest with the 
higher population in round r. Let nL(i,j,r) be the nest with the lower population. For simplicity 
of notation let pn{i, jx) = Pi n H(i, jx)i r ) and pi,(h3-> r ) = p( n H(i, jx), r )■ Define: 




PH(i,j,r) _ 
PL(i,j,r ) 


That is, e(i,j,r) is the relative population gap between the larger and smaller of nests n* and 
nj in round r. 

2 Here it is not important to which nests ants are recruiting and bringing other ants; we just want to reason about 
the value of X'p. Therefore, it is not alarming that ant at is located at nest rij and recruiting with probability p(i, r). 
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Lemma 5.4. For any two nests ni and nj, E [e(i,j, 1)] > 1/(3 (n — 1)). 

Proof. We need to bound the expected difference between two nests after the first round of search¬ 
ing. If nests ni and nj receive the same number of ants after the initial round of searching, then, 
by definition, 1 ) = 0. Let 2 < x < n be the total number of ants that find either nest n* 

or nest nj in the first round of searching. Note that if x = 1 , it is not possible for both nests to 
receive the same number of ants, and if x = 0, then the nests are already empty. The probability 
that both nests receive exactly x/2 ants each is: 


Pr[e(i,j,l) = 0] 


■ x \ ny /2 ny /2 ^ / e2 x \ 1 ^ e < e 2 
x/2) \2y \2/ - fTTy/x) 2 X ~ ir^/x ~ vr\/2 < 3’ 


where in the second step we use the Stirling approximation y/ r 2iTx x+1 ^‘ 2 e~ x < x\ < ex x+1 ^ 2 e~ x . Con¬ 
ditioning on e(i,j, 1) ^ 0, the smallest possible value of e(i,j, 1) is l/(n—1), and so E[e(i, j, l)|e(z,j, 1) / 
0] > l/(n — 1). Therefore, 

E[e(f, j, 1)] > E1)| e(i,j, 1) / 0] • Pr1) / 0] > 7^—[)' □ 


Since e(i,j,r ) is a ratio of two populations, it is not immediately clear how to compute its 
expected change between rounds. However, using a Taylor series expansion we are able to linearize 
this ratio, and then use Lemma 5.3 to compute its expected change. We first show how to use the 
Taylor series expansion in our setting. 


Lemma 5.5. Consider positive numbers xq, x\, yo, and y\. Let Ax = x\ — xq, Ay = yi — yo, and 


A 



Xl_ 

yi 


2a. Then: 
yo 


Ax xqA y 

yo 


vl 


E 

i =0 


-Aj/ 

yo 


(i) 


Proof. We can expand A ^2^ using a Taylor series for f(x,y) = 2. To save notation, we write 
dJdyk-i to denote gj^k-i {x 0 ,yo)- We have: 


A 


df df 1 
yj = fa Ax+ d^ Ay+ 2! 


dx 2 


r ) 2 f 

Ax 2 + 2^xLAxA y + AA Ay 2 


dxdy 


d 2 f 

dy 2 


+ 


3! 


——— Ax^ T ' 

dx 3 


+ ... 


We will show by induction that, for all k: 


1 

~k\ 


d*l 

dx k 


d k f 


Ax k + k n ,, Ax k ~ 1 Ay + ... + —AAy k 


Qkf 


dx k ~ 1 dy 


dy k 


Ax xqA y 

yo 


vl 


—A y 
Vo 


k -1 


( 2 ) 


Plugging into the Taylor series expansion, this immediately gives (1). For k = 1, we differentiate 
f(x,y) = 2 to get our base case: 

df df Ax xqA y 

-frAx + -rf-A y =-— • 

dx dy y 0 y& 


Now for k > 1, by the inductive assumption we have: 


1 


d k ~if 

(k — 1)! [ dx k ~ l 


d k ~i f 

+ ... + ^Ay* 


Ax k ~ 1 1 1 17 J A - fc “ 1 


dy k ~ 


Ax xoAy 
yo vl . 


-A y 

yo 


k 2 Ax(—Ay) k 2 xq(—A y) k 1 


vt 1 


yo 
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Differentiating gives: 


1 

k\ 


^Ax k + ... + ^A y k 
OX K 


d k f 
dy k ‘ 


(-Ay) k 1 (k - l)Ax(-Ay) k 2 ^ kx 0 (-Ay) k 1 ^ 


Vo 


Vo 


Vo +1 


Ax • (—A y) k 1 

xqAx(—A y) k 


Ax 

x 0 Ay~ 

'-Ay' 

O 

5* 

1 

+ v k+1 

Vo 


_ Vo 

Vo . 

. Vo 


k -1 


□ 

For simplicity of notation, write pu = p(nH(i, j,r),r), pl = p(nL(i, j,r),r). Additionally, 
write Ap H = p{n H (i, j, r), r + 2) - p(n H (i, j,r),r), Ap L = p(n L (i, j, r), r + 2) - p(n L (i, j,r),r) and 

A (?r) = p ( (n"( ( S,r},r^2) “ p(""(S,’r),rj ' We CaH 110W a PP^ Lemma 5 ‘ 5 t0 sllOW: 

Lemma 5.6. Let n* and nj be two nests with Api < pj J for some round r € R\. Then, A (jY) P 

1 / a PH _ phAj>l\ 

2 Pl p'i ) ' 

Proof. Note that we are not considering expectations for now, just the actual values. Applying 
Lemma 5.5, we can write A i n terms of A pH and A pl- We have: 


A (= ( ApH - PhApl \ (l - ( ApL \ + f ApL Y _ ( A Pl \ 3 \ 

\Pl ) V PL P 2 L J \ \ Pl J V Pl J \ Pl ) 

The value j s strictly less than 1 by the assumption that A pl < Pl- So: 


Next, we use the value of A (ph/pl) in order to calculate the expected change in e(i,j,r ) for 
two nests rij and nj after one round of recruiting. 

Lemma 5.7. Let n{ and nj be two nests with A pl < Pl, p{i,i") > l/(dk) and p(j,r) > 1 /(dk) for 
some round r G R\. Then, E [e{i,j,r + 2)] > (1 + 1/(2 dk)) E [e{i,j,r)\. 

Proof. First we show that: 



A plY _ l_ 

Pl J 1 + a pl/pl “1 + 1 


> 


= 1/2. 


□ 



e{i,j,r+ 2) > 


p(n H (i,j,r),r + 2) 
p(n L (i, j,r),r + 2) 


- 1. 


This is true because, if the larger of the two nests in round r remains larger in round r + 2, then 
nuif, j, f) = nn(i,j, r + 2) and nL{i,j , r) = nL(i,j, r + 2), so, by the dehnition of e, e(i,j, r + 2) = 


P(n-H(i,j , r )> r + 2 ) 

p(nL(i,j,r),r+2) 


— 1. On the other hand, if the two nests flip positions, then 


P(n H (i,j,r),r+ 2) 
p(nL(i,j,r),r+ 2) 


— 1 must 


be smaller than e{i,j,r + 2) since it contains the population of the larger nest in the denominator 
of the fraction. Fixing p(i,r ) and p(J,r), we have: 


E[e(z,j,r + 2)] > E 


p(n H (i,j,r),r + 2) 
p{n L {i,j,r),r + 2) 



e(i,j,r)+E 
e(i,j,r) + E 


'p{n H (i,j,r),r + 2) ,. . ' 

—-;----— 1 — e[i, 7, r) 

.p( n L(hj,r),r + 2) 

' p{n H (i,j,r),r + 2) _ p(n H {i, j,r),r) 

_p(n L {i, j,r),r + 2) p(n L (i,j,r),r) 
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So we have reduced our question of how we expect the absolute difference between the nests to 
change in a single round to the question of how we expect the ratio of population to change. 

By Lemma 5.6, it follows: 

e(i,j,r + 2)\ > e(i,j,r) + ^ fE 

By Lemma 5.3 we have: 

p L E[Ap H ] -p H E[Ap L \ = PLP 2 H fi{n H ,r) - PLPH&{riH,r) ~ PHP 2 L fi(n L ,r) + PLPHf2{ji L ,r) 

= PlPh ( fi H ,r) - &{nH,r) - pUi ( n L , r) + £ 2 ( n L , r )] 

> PlPh£,i(ji l , r ) • ( p H - p L ) > £,plPh{ph ~ Pl), 


A PH PH^PL 


PL 


pi 


> 4 iJ, r) + i ( WE[Am]-mE[A p L ] 

2 V Pl 


where the inequality follows from the guarantee of Lemma 5.2 that PL£i( n L> r ) — £ 2 (ut,,r) < 

ig that e(i,j , r) = 1 

1 (PlPh (ph ~ Pl) 


PLfi{nH,r) — £2 {riH,r). Noting that e(i,j,r ) = Ph Pi Pl we have 


E[e(i,j,r+ 2)] > e(i,j,r) + 


2 V 


v\ 


-)>e(i,J.r)(l + ^)> £ (M,r)(l + jL). 


□ 


5.2.3 Changes in Populations of Nests over 0(k log n) Rounds 

We now show that once the population of a nest is very small, it will quickly drop to zero. 


Lemma 5.8. Let ni be any nest with p(i, r ) < 1 /(dk) in some round r € R\. Then, with probability 
at least 1 — l/n c+4 , p(i,r') < 1 /(dk) for all rounds r' € [r,r + 64(c + 4)/clogn)]. 

Proof. It suffices to show that if p(i, r) < 1 /(dk) thenp(i, r+2) < 1 /(dk) with probability 1—l/n c+5 . 
We can then union bound over all 0(klogn) rounds to get that, with probability 1 — l/n c+4 , 
p{i,r') < 1 /(dk) for all rounds r' € [r, r + 64(c + 4)k log n)]. 

Consider two possible cases: 


Case 1: p(i,r) < 1/(2 dk). Then, even if all ants successfully recruit, the nest cannot more than 
double in size, so p(i, r + 2) < 1 /(dk). 

Case 2: p(i,r) > 1/(2 dk). The expected number of ants to recruit to nest n, is n ■ p(i,r) 2 > 
n/(4d 2 k 2 ). By the assumption that k < ^/n/(8d 2 (c + 6) log n) = 0{y/n/ log n), this value is 
in fl(logn). By a Chernoff bound, with probability at least 1 — l/n c+6 , at most 2 n ■p(i , r) 2 = 
n/{2d 2 k 2 ) ants choose to recruit. 

Similarly, the expected the number of ants to be recruited from nest rij is: 


np(i, r) (1 — p(i, r)) L 2 (r) > n 


1 

dk 


1 - 


2 dk 


1 n 

k-dif 2 


1 - 


2 dk 


where in the second step we use the fact that 1 /(dk) > p(i,r) > 1/(2 dk) and X 2 (r) > 1 /k. 
This value is also in fl(logn). 

Next, we would like to apply a Chernoff bound to the number of recruited ants from nest 
np, however, recruitments are not independent, so we need to do a simple domination trick, 
similarly to the proof of the lower bound. 

Let random variable Yf 1 have value 1 with probability 1/A;(1 — 1/(2 dk)), and let Y r be the 
sum of n/(dk) such independent random variables. So, E[y r ] = (n/(dk 2 ))( 1 — 1/(2 dk)), and 
by a Chernoff bound, P[Y r < (3/2)(n/(d/c 2 )(l — 1/(2 dk))\ < n^( c+6 \ We can see that: 


17 











P[a is recruited | the choices of all other ants] > (1 ~p{i,r))E 2 (r) >11 — 


1 \ 1 


2 dk J k 


= P[Y r a = 1] 


Therefore, by Lemma A.l, P[ff of ants recruited from nest n* < x] < P[Y r < x\ for some 
x < n/(dk). Plugging in x = (3/2)(n/(dfc 2 )(l — 1/(2 dk)), it follows that: 

P[ff of ants recruited from nest ra* < (3/2)(n/(d/c 2 )(l — 1/(2 dk))\ < n~( c+6 \ 

Therefore, with probability at least 1 — l/n c+6 , at least (3/2)(n/(dk 2 )(l — 1/(2 dk)) ants will 
be recruited from nest rii. Recall that, with probability at least 1 — l/n c+6 , at most n/{2d 2 k 2 ) 
ants from nest n* choose to recruit. Since d > 64, it is easy to see that the number of ants 
recruiting is smaller than the number of ants recruited away from nest n*, so with probability 
at least 1 — l/n c+ , the total population of rtj decreases, indicating p(i,r + 2) < 1 /(dk). 

So overall, if p(i,r ) < l/(dk), then, with probability at least 1 — l/n c+4 , p(i,r') < 1 /(dk) for 
all r' £ [r,r + 0(k log n)\. □ 

Lemma 5.9. Let rii be a nest with p(i,r ) < 1 /(dk) in some round r € R\. Then, with probability 
at least 1 — l/n c+3 , c(i,r') = 0 for r' = 64(c + A)k log n). 

Proof. We calculate the expected change in the number of ants in nest rii by first calculating E[X“] 
for some ant a in nest n* and some round r € R\ such that p(i,r) < 1 /(dk). 

Suppose ant a does not recruit. Let random variable Yf 1 have value 1 if ant a' successfully 
recruits ant a, and value 0 otherwise. By Lemma 2.1, it follows that E[T r a ] > p(if, r) / (16n) where 
i! = £(a',r). Therefore: 


E 


E F “' 


a'^a 


E E [ y “'] 


a'^a 


^itE E *'.’■) 

j=1 {a'jt=a\£(a' ,r)=j} 


2 E 

\i=l {a'\£(a',r)=j} J 

a ^ (E"p(j.o s j 


dk ~ 16 dk 


Therefore, the expected value of Xf is: 


, WI £ 5-( 1 -S i ‘ 


E E 

{a’\i{a’,r)^i} 


= L_L_L1 (EE _ E) < _2_, 

dk V dk J l 16 dk J 64 k 


where the last step follows from the assumption that d > 64. Therefore, 

E\p(i,r + 2)]=p(i,r) + ^ E E [ X “] < p(i, r)(l + E[X“]) < p(i, r) ^1 - 


{a|£(a,r)=i} 
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By Lemma 5.8, p(i,r' ) < 1 /(dk) for all r' G [r, r + 64(c + 4)Llogn)]. Therefore, for r / = 
64(c + 4)A:logn), it is true that E [p(i,r')\ < l/n c+3 and, by a Markov bound, nest n* has at 
least one ant with probability at most l/n c+3 . Therefore, with probability at least 1 — l/n c+3 , 
c(i, r') = 0. □ 

Next, we show that, for any pair of nests both with population proportions bounded below by 
1 /(dk), we can use Lemma 5.7 to argue that the populations of these nests diverge quickly. As 
soon as a nest drops below the 1 /(dk) threshold we can use Lemma 5.9 to show that it will not be 
the winning nest. 

Lemma 5.10. Letrii andrij be two nests with q(i) = q(j) = 1. Then, forr' = (6d+64(c+6))L logn, 
with probability at least 1 — l/n c+2 , at least one of the following is true: c(i,r') = 0 or c(j,r') = 0. 

Proof. Note that if at any point a nest has no ants in it, it remains having no ants thereafter. We 
consider two possible cases based on how many ants are in each nest in each round r € [1, 6 dk logn)]: 

Case 1: In some round r G [1,6 dk logn)] either n* of nj has fewer than n/(dk) ants. Then, by 
Lemma 5.9, with probability at least 1 — l/n c+3 , this nest has no ants committed to it after 
64(c + 4 )k log n) rounds. 

Case 2: In every round r G [1,6 dk log n)] both n* and rij have at least n/(dk) ants. We are going 
to show that, with probability at least 1 — l/n c+3 this case does not happen. Suppose, in 
contradiction, that in every round r € [1,6 dk log n)] both n* and rij have at least n/(dk ) ants. 
First, we want to show that Api < pl, so that we can apply Lemma 5.7. Since both nests 
have at least n/ (dk) ants, it follows that pl > 1 /(dk), so the number of ants recruiting for that 
nest is at least n/(d 2 k 2 ) = fl(logn), by our assumption that k < yjn/ (8 d 2 (c + 6) log n) = 
0(y/n/ logn). By a Chernoff bound, with probability at least 1 — l/n c+3 , the number of 
recruiting ants is at most 2 n/(d 2 k 2 ). This indicates that, with probability at least 1 — l/n c+3 , 
not all ants from the nest with lower population recruit, and so A pl < Pl- 

This means that by Lemma 5.7, 

E[e(i,j,dk log n)}> ( 1 + E[e(i, j, 1)] > n 3 • E[e(z, j, 1)]. 

By Lemma 5.4, E[e(*,j, 1)] > l/(3(n — 1)). Therefore, after 6 dklogn rounds we have 
E[e(i, j, 6dk logn)] > n. However, this is a contradiction to the fact that e(i,j, 6dk log n) < 
n — 1, by definition, and consequently E[e(?', j , 6 dk log n)] < n — 1. 

Each case holds with probability at least 1 — l/n c+3 , so union bounding them, we get that, with 
probability at least 1 — l/n c+2 , either c(i,r') = 0 or c(j,r') = 0. □ 

Theorem 5.11. With probability at least 1 —l/n c , Algorithm 3 solves the HOUSEHUNTING problem 
in O (k logn) rounds. 

Proof. By our bound on k (k < i/ n/(8d 2 (c + 6) logn) = 0(\Jn/ logn)), with probability at least 
1 — l/n c+1 , in the first round of the algorithm (a round of searching), at least some ant will arrive 
at a nest with quality 1. Further, since only ants at nests with quality 1 choose to recruit, at any 
point, at least one ant is committed to a good nest. 

By Lemma 5.10, with probability at least 1 — l/n c+2 , for each pair of nests rij and rij , at 
least one nest contains no ants by the end of 0(k\ogn) rounds. Union bounding over all, at most 
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( 2 ) < k 2 < n (by the bound on k), pairs of good nests, we conclude that, with probability at least 
1 — l/n c+1 , after 0(k logn) rounds, each pair contains at least one nest with no ants. This can 
only be true if all nests have no ants (not possible) or if all ants are located at one good nest. 

Finally, union bounding the initial search phase and the subsequent competition between nests, 
we get that with probability at least 1 — l/n c , the house hunting problem is solved in O(k logn) 
rounds. □ 

6 Discussion and Future Work 

Extensions to the Model For the sake of analysis, we have made many simplifying assumptions 
about the house-hunting process. We are confident that many of these assumptions can be weakened 
to make the model more realistic and natural. Some obvious modifications include assuming ants 
know only an approximation of n, allowing values of k larger than 0(n/logn), and allowing non¬ 
binary nest qualities, variability in the ants’ quality sensing, along with some measure of algorithmic 
performance based on the quality of the chosen nest. Distinguishing between direct transport and 
tandem runs may also be interesting, paired with a more fine-grained runtime analysis. 

Additionally, real ants can only assess nest quality and population approximately. For example, 
they may estimate nest size (one measure of quality) by randomly walking within the nest and 
counting how many times they cross their previous path [20]. They seem to estimate nest population 
by measuring encounter rates with other ants, with a higher encounter rate indicated a higher 
population at the nest [14, 21]. Adding noisy measurements to our model and designing algorithms 
that handle this noise would be a very interesting future direction. It may even be possible to 
explicitly model lower level behavior and implement subroutines for nest assessment, population 
measurement, recruitment, and search which give various runtime and error guarantees. 

Extensions to the Algorithms We believe that Algorithm 3 may be a good starting point for 
work on more realistic house-hunting models. Below we discuss some interesting possible extensions 
to the algorithm. Some seem to simply require a more involved analysis, while others seem to require 
trade-offs in the algorithm’s running time and its level of simplicity. 

Improved running time: The 0{k\ogn) runtime of Algorithm 3 is required because, on 
average each nest initially contains n/k ants, so ants only recruit with probability 1 /k. 0(k) time 
is required to amplify population gaps by a constant factor. Ideally, ants would all recruit with 
a probability lower bounded by a constant, but still linearly dependent on the nest populations. 
This would allow convergence in 0(logn) rounds. If ants keep track of the round number, they 
can map this to an estimate k[r ) of how many competing nests remain, allowing them to recruit at 
rate 0(c(i,r)/n ■ k(r)). We believe that such a strategy should yield a relatively natural algorithm 
converging in 0(log c n) rounds. 

Non-binary nest qualities: Assuming a real-valued nest quality in the range (0,1) affects 
the correctness of Algorithm 3 because ants no longer have the notion of a good nest. However, it 
should be possible to incorporate the quality of the nest into the recruitment probability in order 
make the algorithm converge to a high-quality nest, without significantly effecting runtime. 

Approximate counting, nest assessment, and knowledge of n: Since the analysis of 
Algorithm 3 does not require each ant to recruit with a specific probability, but rather that the 
total number of ants choosing to recruit from a nest is proportional to its population, it should be 
resilient to noisy quality and population measurements. As long as ants have unbiased estimators 
of these values, we believe that Algorithm 3 should remain correct, perhaps with some runtime cost 
dependent on estimator variance. 
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Fault tolerance: Similarly, Algorithm 3 should support some degree of fault tolerance. A 
small number of ants suffering from crash-faults or even malicious faults, should not affect the 
overall populations of recruiting ants and the algorithm’s performance. 

Asynchrony: Finally, note that Algorithm 3 currently works in synchronous rounds and relies 
on that assumption to get the correct number of ants at a given nest. However, we believe that, as 
long as the distribution of ants in candidate nests throughout time stays close to the distribution 
in the synchronous model, Algorithm 3 can be extended to work in a partially-synchronous model, 
potentially at the cost of some extra running time. 


References 

[1] Karl Johan Astrom and Richard M Murray. Feedback systems: an introduction for scientists 
and engineers. Princeton university press, 2010. 

[2] Fan Chung, Shirin Handjani, and Doug Jungreis. Generalizations of Polya’s urn problem. 
Annals of combinatorics, 7(2): 141—153, 2003. 

[3] Alejandro Cornejo, Anna Dornhaus, Nancy Lynch, and Radhika Nagpal. Task allocation in 
ant colonies. In Distributed Computing, pages 46-60. Springer, 2014. 

[4] Adam L Cronin. Synergy between pheromone trails and quorum thresholds underlies consensus 
decisions in the ant myrmecina nipponica. Behavioral Ecology and Sociobiology, 67(10): 1643— 
1651, 2013. 

[5] Benjamin Doerr and Anne Auger. Theory of Randomized Search Heuristics. World Scientific, 

2011. 

[6] A Dornhaus and NR Franks. Colony size affects collective decision-making in the ant tem- 
nothorax albipennis. Insectes sociaux, 53(4):420-427, 2006. 

[7] Yuval Emek, Tobias Langner, Jara Uitto, and Roger Wattenhofer. Solving the ANTS problem 
with asynchronous finite state machines. In Proceedings of the fist International Colloquium 
on Automata, Languages, and Programming (ICALP), 2014. 

[8] Ofer Feinerman, Bernhard Haeupler, and Amos Korman. Efficient and noise-resilient rumor¬ 
spreading in large, anonymous populations. arXiv preprint arXiv:1311.3f25 , 2013. 

[9] Ofer Feinerman and Amos Korman. Memory lower bounds for randomized collaborative search 
and implications for biology. In Distributed Computing, pages 61-75. Springer, 2012. 

[10] Ofer Feinerman and Amos Korman. Theoretical distributed computing meets biology: A 
review. In Distributed Computing and Internet Technology, pages 1-18. Springer, 2013. 

[11] Ofer Feinerman, Amos Korman, Zvi Lotker, and Jean-Sebastien Sereni. Collaborative search on 
the plane without communication. In Proceedings of the 2012 ACM symposium on Principles 
of distributed computing, pages 77-86. ACM, 2012. 

[12] Michael J Fischer, Nancy A Lynch, and Michael S Paterson. Impossibility of distributed 
consensus with one faulty process. Journal of the ACM (JACM), 32(2):374—382, 1985. 


21 


[13] Nigel R. Franks, Anna Dornhaus, Charlotte S. Best, and Elizabeth L. Jones. Decision making 
by small and large house-hunting ant colonies: one size fits all. Animal Behaviour, 72(3):611 
- 616, 2006. 

[14] Deborah M Gordon. Ant encounters: interaction networks and colony behavior. Princeton 
University Press, 2010. 

[15] Christiane IM Healey and Stephen C Pratt. The effect of prior experience on nest site evalu¬ 
ation by the ant Temnothorax curvispinosus. Animal Behaviour, 76(3):893-899, 2008. 

[16] Richard Karp, Christian Schindelhauer, Scott Shenker, and Berthold Vocking. Randomized 
rumor spreading. In Foundations of Computer Science, 2000. Proceedings. 41st Annual Sym¬ 
posium on, pages 565-574. IEEE, 2000. 

[17] Leslie Lamport. The part-time parliament. ACM Transactions on Computer Systems (TOCS), 
16(2):133-169, 1998. 

[18] Christoph Lenzen, Nancy Lynch, Calvin Newport, and Tsvetomira Radeva. Trade-offs between 
selection complexity and performance when searching the plane without communication. In 
Proceedings of the 2014 ACM symposium on Principles of distributed computing, 2014. 

[19] E Mallon, S Pratt, and N Franks. Individual and collective decision-making during nest site 
selection by the ant Leptothorax albipennis. Behavioral Ecology and Sociobiology, 50(4):352- 
359, 2001. 

[20] Eamonn B Mallon and Nigel R Franks. Ants estimate area using Buffon’s needle. Proceedings 
of the Royal Society of London. Series B: Biological Sciences, 267(1445):765-770, 2000. 

[21] S.C. Pratt. Nest site choice in social insects. In Michael D. Breed and Janice Moore, editors, 
Encyclopedia of Animal Behavior, pages 534 - 540. Academic Press, Oxford, 2010. 

[22] Stephen C Pratt. Quorum sensing by encounter rates in the ant Temnothorax albipennis. 
Behavioral Ecology, 16(2):488-496, 2005. 

[23] Stephen C Pratt, Eamonn B Mallon, David J Sumpter, and Nigel R Franks. Quorum sensing, 
recruitment, and collective decision-making during colony emigration by the ant Leptothorax 
albipennis. Behavioral Ecology and Sociobiology, 52(2):117-127, 2002. 

[24] Stephen C Pratt and David JT Sumpter. A tunable algorithm for collective decision-making. 
Proceedings of the National Academy of Sciences, 103(43):15906-15910, 2006. 

[25] Takao Sasaki and Stephen C Pratt. Emergence of group rationality from irrational individuals. 
Behavioral Ecology, 22(2):276-281, 2011. 

[26] Takao Sasaki and Stephen C Pratt. Ants learn to rely on more informative attributes during 
decision-making. Biology letters, 9(6):20130667, 2013. 

A Math Preliminaries 

Lemma A.l. Let X\, - ■ ■ ,X n be arbitrary binary random variables. Let X(,--- ,X* be random 

variables that are mutually independent and such that for all i, X* is independent of X i, • • • , . 

Assume that for all i and all x i, • • • ,Xi- 1 € {0,1}, 
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P[Xi = l|Xi = x ir • ■ = Si_!] > P[X* = 1] 


Then, for all k > 0, we have, 


p 

n 

J2 X i<k 

< p 

n 

J2 x i <k 


_i =1 


_i=1 


and the latter term can be bounded by Chernoff bounds for independent random variables. 
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